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The matter density distribution inside a hadron is evaluated using gauge-invariant correlation functions within 
the quenched and the unquenched theory. Comparison with the charge density distribution suggests that hadron 
deformation is a consequence of the relativistic motion of the quarks. 

1. Introduction 



The charge and matter density distributions 
encode detailed information on hadron structure. 
In the non-relativistic limit both charge and mat- 
ter density distributions reduce to the square 
of the wave function. Therefore differences be- 
tween them give a measure of relativistic effects. 
Hadron deformation can be searched via the mat- 
ter density distribution as well as via the charge 
density distribution. There are strong experimen- 
tal indications from measurements of the electro- 
magnetic transition form factors in 7* A*" to A [1] 
that the nucleon is deformed. It is thus interest- 
ing to address the issue of hadron deformation 
within lattice QCD and understand its physical 
origin. Comparison of matter and charge density 
distributions can help discriminate among possi- 
ble mechanisms. A number of phenomenological 
studies use the bag-model to describe hadron de- 
formation. For this reason we compare the charge 
and matter distributions to the bag-model results. 

2. Observables 

We evaluate equal time two-current correlators 

Cr(r,t) = J {h\ff{r' + r,t)j^(r\t)\h) (1) 

with the current operator given by the normal 
order product jp(r, —: u{r,t)Tu{r,t) : and 
r — 7o, 1 for the charge and matter correlators 
respectively. We note that no gauge ambiguity 
arises in this determination of hadron wave func- 
tions unlike Bethe-Salpeter amplitudes. In the 
case of baryons two relative distances are involved 
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and three current insertions are required. How- 
ever we may consider integrating over one relative 
distance to obtain the one-particle density distri- 
bution that involves two current insertions and it 
is thus evaluated via Eq. 1. 

3. Results 

All the results presented here have been ob- 
tained on lattices of size 16"^ x 32. We ana- 
lyze, for the quenched case, 220 NERSC config- 
urations at P — 6.0 for pion to rho mass ratio 
m^/mp = 0.88, 0.84, 0.78 and 0.70 and, for the 
unquenched, 150 SESAM configurations [4] for 
K = 0.156 and 200 for k ~ 0.157 with niT^/nip = 
0.83 and 0.76 respectively. The physical volume 
of the quenched and unquenched lattices is ap- 
proximately the same. 

In Eq. 1 the current couples to the quark at 
fixed time separation t from the source, which 
must be large enough to sufficiently isolate the 
hadron ground state. Using local and Wuppertal 
smeared sources we first checked that we obtain 
the same value for the effective mass at i = 8a. 
This value agrees with the one extracted using 
data with Dirichlet boundary conditions utilizing 
the whole time extent of our lattice. In contrast, 
wall sources even at maximal separation t = 8a 
still show a sizable excited states contamination. 
In addition, we checked, by comparing data at dif- 
ferent t, that higher momenta are sufficiently sup- 
pressed in the current-current correlators since 
summing over the spatial volume of the source is 
not possible. Also, results obtained using Wup- 
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Figure 1. C/(r) for different quark masses. The 
errors bars are not shown for clarity. 

pertal smearing which suppresses high momenta 
and local sources which allow them, are in agree- 
ment at i = 8a. This means that contamination 
from high momentum states is negligible at f = 8a 
allowing us to analyze standard full QCD configu- 
rations that employ anti - periodic b.c. Since, for 
our parameters, local sources produce the same 
results for Cr{r) as smeared ones and carry no 
gauge noise they are preferable for this study. 

The quenched matter density distribution 
shown in Fig. 1 broadens for the pion as the 
quarks become lighter. For the rho, the nucleon 
and the A"*" essentially no variation is seen over 
the range of naive quark masses ~ 300—100 MeV. 

1.00 





1 1 1 

charge 


II 

\ matter 






/c = 0.154 


V /c = 0.154 




\z - 








\z 









r/a r/a 

Figure 2. Cr(0,0,2;) (upper Z-curve) and 
Cr(x,0,0) and Cr(0,y,0) (lower X,y-curves). 

In a previous study of the charge density dis- 
tribution [2,3] we obtained a non-zero quadrupole 
moment for the rho or a 2 — x asymmetry where 
the 2-axis is taken along the spin of the rho as 
shown in Fig. 2. An angular decomposition of 
the wave function, as shown in Fig. 3, corrobo- 
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Figure 3. Decomposition into L = and L = 2 
contributions. 

rates a non-zero charge deformation by producing 
a non-zero L = 2 component. No such deforma- 
tion is seen for the matter density distribution in 
Figs. 2,3. Since the matter and charge operators 
have the same non-relativistic limit, this strongly 
suggests that hadron charge deformation is a rela- 
tivistic effect. This result has strong implications 
for the validity of various phenomenological mod- 
els used in the study of nucleon deformation and 
•j^N — > A form factors. 
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Figure 4. (a) Quenched and (b) unquenched mat- 
ter density distributions. 

Fig. 4 compares quenched and unquenched re- 
sults for C/(r) at TOt/top ~ 0.83. Unquenching 
increases C/(r) at short distances in the case of 
the pion and the rho whereas for the baryons no 
significant changes are seen. Both quenched and 
unquenched matter distributions show a faster 
fall off than the charge distribution. Further- 
more, whereas unquenching tends to lead to a 
small increase in the rho charge asymmetry, it 
has no effect on the matter density distribution. 
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This again suggests a relativistic origin for the 
rho charge deformation. 

4. Comparison of lattice and bag model re- 
sults 

We consider only the lowest mode of the free 
Dirac field in a spherical bag of radius R that is 
chosen so as to minimize the mass of the hadron 
under consideration. For the comparison pre- 
sented here, we use the lattice values for the 
masses of the rho, the nucleon and the A to fix the 
bag parameters, B,Zo and atag, using as an in- 
put the naive quark mass. Expanding the quark 
fields in terms of bag eigenmodes we obtain for 
the charge and matter density distributions 

(/(r)2+.g(r)2) 
(/(r')2+.9(r')2) 

{f{r'?~9{r'?) 



Ji{r)jf{r' 



= C' 



(2) 



(3) 



with C, C constants. M denotes mesons and B 
baryons. These expressions show explicitly that 
the difference between charge and matter density 
distributions is due to the opposite sign of the 
lower components of the Dirac spinor i.e. a rela- 
tivistic effect. 
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Figure 5. Quenched charge and matter distribu- 
tions at K = 0.153 and bag model results (dashed 
line). The solid line is a fit to the form exp(— mr). 

From Fig. 5 we see that the bag model fails 
to reproduce the correct radial dependence of the 
lattice results for the charge and matter density 



distributions and that the Ansatz exp(— mr) pro- 
vides the best description of the data. A measure 
of the width of the distributions is provided by 
the root mean square (rms) radius. The ratio of 
the charge to the matter rms radius in the bag- 
model at K = 0.153 and 0.154 is: 1.16 and 1.17 
for the rho and 1.16 for the nucleon. The cor- 
responding ratios for the lattice data are 1.15(3) 
and 1.16(8) for the rho and 1.18(4) and 1.20(8) for 
the nucleon. Thus, despite the failure of the bag 
model in describing the individual radial shape 
of the distributions it produces reasonable results 
for the relative widths of the charge to the matter 
distribution. Both lattice and bag model consis- 
tently predict a broader charge than matter dis- 
tribution, with a very weak mass dependence. 

5. Conclusions 

For the lattice parameters used in this study, 
local sources are suitable interpolating fields for 
the evaluation of density distributions, since they 
have less gauge noise than smeared ones and the 
temporal extent of our lattice is large enough 
to obtain the ground state. Quenched and un- 
quenched results for both the charge and mat- 
ter density distributions show insignificant dif- 
ferences. The charge density distribution is, in 
all cases, broader than the matter density. For 
baryons, the lattice indicates a charge radius 
about 20% larger than the matter radius. This 
effect is well reproduced by the bag model. The 
deformation seen in the rho charge distribution 
is absent in the matter distribution, both in the 
quenched and the unquenched theory. This ob- 
servation suggests a relativistic origin for the de- 
formation. 
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